EXISTENCE OF 3D STRONG SOLUTIONS FOR A SYSTEM 
MODELING A DEFORMABLE SOLID IN A VISCOUS 
INCOMPRESSIBLE FLUID* 
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Abstract. In this paper wc study a coupled system modeling the movement of a deformable solid 
immersed in a fluid. For the solid we consider a given deformation that has to obey several physical 
constraints. The motion of the fluid is modeled by the incompressible Navier-Stokes equations in a 
time-dependent bounded domain of R 3 , and the solid satisfies the Newton's laws. Our contribution 
consists in adapting and completing some results of [16] in dimension 3, in a framework where the 
regularity of the deformation of the solid is limited. We rewrite the main system in domains which 
do not depend on time, by using a new means of defining a change of variables, and a suitable change 
of unknowns. We study the corresponding linearized system before setting a local-in-time existence 
result. Global existence is obtained for small data, and in particular for deformations of the solid 
which are arbitrarily close to the identity. 
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1. Introduction. 

In this paper we arc interested in a deformable solid immersed in a viscous incom- 
pressible fluid in dimension 3. The domain occupied by the solid at time t is denoted 
by S(t). We assume that S(t) C O, where O is a bounded regular domain. The fluid 
surrounding the solid occupies the domain O \ S{t) = F(t). 



O = T{t) U S{t) C 




1.1. Presentation of the model. 

The movement of the solid in the inertial frame of reference is described through the 
time by a Lagrangian mapping X$, so 

S(t) = X s (S(0),t), t>0. 
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The mapping Xg(-,t) can be decomposed as follows 

X s (y,t) = h(t)+H(t)X*(y,t), y € 5(0), 

where the vector h(t) describes the position of the center of mass and R(i) is the 
rotation associated with the angular velocity of the solid. More precisely, u> and R 
are related to each other through the following Cauchy problem 

(0 -W3 uj 2 
Ll>3 — U)\ 

—U)2 Ml 

The couple (/i(i),R(i)) describes the position of the solid and is unknown, whereas 
the mapping X*(-,t) can be imposed. This latter represents the deformation of the 
solid in its own frame of reference and will constitute the main datum of the problem. 
When this Lagrangian mapping X*(-,t) is invertible, we can link to it an Eulerian 
velocity w* through the following Cauchy problem 

dX* 

— (y,t)=w*(X*(y,t),t), X*(y,0) = y-h(0), y G 5(0). 

If Y*(-,t) denotes the inverse of X*(-,t), we have 

dX * 

w*(x*,t) = — (Y*{x*,t),t), x* e S*(t) = X*(S{0),t). 

This Eulerian velocity w* can also be considered as a datum defining the way the solid 
is deforming itself. Considering X* - or w* - as a datum is equivalent to assuming 
that the solid is strong enough to impose its own shape. 

The fluid flow is described by its velocity u and its pressure p. For w* satisfying a set 
of hypotheses given further, we aim at proving the existence of strong solutions for 
the following coupled system 

On 

— -vAu+(u-V)u + Vp = 0, ieJ(i), fe(0,T), (1.1) 
ot 

div u = 0, i£ T(t), t e (0, T), (1.2) 

u = 0, xedO, te(0,T), (1.3) 



u 



h'(t)+(j{t)A(x-h(t))+w{x,t), xe dS(t), te(0,T), (1.4) 

Mh"(t) = - [ <j{u,p)ndT, te{0,T), (1.5) 

(Iu>)'{t) = -[ (x- h(tj) Aa{u,p)ndT, t e (0,T), (1.6) 

JdS(t) 



u{y, 0) = u {y), y e J"(0), h(0) = h e M 3 , h'{0) =h x e M 3 , uj(0) = uj G M 3 , 

(1.7) 

where 



S(t) = h{t)+R(t)X*(S(0),t), J r (t)=0\S(t), (1.8) 
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and where the velocity w is defined by the following change of frame 



(x,t)=B.(t)w*(R(ty(x-h(t)),t), i6 5(i). (1.9) 



Without loss of generality, we can assume that h = 0, for a sake of simplicity. The 
symbol A denotes the cross product. The linear map uj A • can be represented by 
the matrix S(uj). In equations (1.5) and (1.6), the mass of the solid M is constant, 
whereas the moment of inertia tensor depends on time, as 

I(t)= [ p s (x,t) (\x - h(t)\ 2 I R3 - (x - h(t)) ® (x - h{t))) dx. 
Js(t) 

The quantity p$ denotes the density of the solid, and obeys the principle of mass 
conservation 

»<*'<»■'>■'> = 9ES(0) ' 

where VXg is the Jacobian matrix of mapping X$. We can define 

_ Ps or*(x;t),o) 

9 {X > t} - det(VX*(Y*(x*,t),t)Y xe6 W- 
For a sake of simplicity we assume that the solid is homogeneous at time t = 0: 

ps(y,o) = ps>o. 

In system (1.1)— (1.9), v is the kinematic viscosity of the fluid and the normalized 
vector n is the normal at dS(t) exterior to J-{t). It is a coupled system between 
the incompressible Navier-Stokes equations and the Newton's laws. The coupling is 
in particular made in the fluid-structure interface, through the equality of velocities 
(1.4) and through the Cauchy stress tensor given by 

a(u,p) = 2vD(u) -<pld = v (Vu + (Vii) T j - p Id. 

We assume that the deformation X* satisfies a set of hypotheses: 
HI For all t £ [0,T], X*(-,t) is a C°°-diffeomorphism from 5(0) onto S*(t). 
H2 In order to respect the incompressibility condition given by (1.2), the volume of 
the whole solid is preserved through the time. That is equivalent to say that 

f f OX* 

/ w* ■ n&Y = / — — ■ (cofVX*) ndT = 0, (1.10) 

JdS*(t) JdS{0) Ot 

where coiVX* denotes the cofactor matrix of VX* . 
H3 The linear momentum of the solid is preserved through the time, that means 

p*(x*,t)w*(x*,t)dx* = PS / (y,t)dy = 0. (1.11) 

S*(t) Js(p) at 



H4 The angular momentum of the solid is preserved through the time, that means 

f OX* 
p*(x*,t)x*Aw*(x*,t)dx* = p s / X*(y,t) A -——(y,t)dy = 0. 

5*(t) JS(0) at 

(1.12) 
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Imposing constraints (1.11) and (1.12) enables us to get the two following constraints 
on the velocity w 

[ p s (x,t)w(x,t)dy = 0, (1.13) 
J sit) 

/ ps(x, t)(x — h(t)) A w(x, t)dy = 0. (1.14) 
J sit) 

As equations (1.5) and (1.6) are written, equalities (1.13) and (1.14) are already 
assumed in system (1.1)-(1.9). Hypotheses H3 and H4 arc made to guarantee the 
self-propelled nature of the motion of the solid, that means no other help than its own 
deformation enables it to move in the fluid. By the undulatory motion induced by its 
own internal deformation, the solid imposes partially, through w, the nonhomogeneous 
Dirichlet condition (1.4). The latter induces the behavior of the environing fluid 
through (1.1)— (1.3), and thus the response of the fluid - given by a(u,p)n on the 
interface - enables the whole solid to be carried, regarding to the ordinary differential 
equations (1.5) and (1.6). The other part of the interaction consists in the fact that 
domains occupied by the fluid and the solid change through the time, as follows 

S(t) = h(t) + K{t)S* (t), T(t) = 0\ S(t). 

1.2. Main result and contributions. The main result we state in this paper 
is Theorem 6.1, that we give as follows: 

Theorem 1.1. 

Assume that < dist(<S, dO), and that uo £ H 1 (J r ) satisfies 

div mo = in J 7 , uq = on dO, uo(y) = hi + luq A y on dS. 

Assume that X* satisfies the hypotheses HI— H4 and that the displacement X* —Ids is 
small enough in Wq^O, oo\S), for some integer m > 3. Assume also that \\uo\\ni (f) , 
\hi\^3 and |wo|r3 are small enough. Then problem (1.1) —(1.9) admits a unique global 
strong solution (it,p, h' , uj) in 

U{Q,oo\T(t)) x L^oojH 1 ^))) x H^O, oc;]& 3 ) x 1^(0, oo;R 3 ). 



This type of problem has been studied in [16] in 2 dimensions, in the case where 
no limitation was supposed on the regularity of the mapping X* . In particular global 
existence is obtained without smallness assumption on the data. We extend this result 
to dimension 3 in a framework where the regularity of the mapping X* is limited, 
which has not been done yet for this system, as far as we know. The strategy for 
proving the existence of strong solutions is globally the same as the one used in [17, 7] 
(for rigid solids), [16], or even in [2] for instance: We first define a change of variables 
which enables us to set a change of unknowns whose the space domain of definition 
does not depend on time anymore. Then we write the nonlinear system that have to 
satisfy the new unknowns, and we study the linearized system associated with. Then 
a local-in-time existence result is proven by a fixed point method, and the global ex- 
istence is obtained by writing appropriate energy estimates. 
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In addition to the technical difficulties induced by the framework of dimension 3, 
the originality of our approach lies in the fact that we have to develop new means 
of handling a deformation of the solid which is limited in regularity. First, a new 
method is developed in order to define a change of variables in the fluid part. Indeed, 
the method introduced in [17] cannot be applied anymore, or at least not so straight- 
forwardly anymore; The way proceed is more direct and more adapted for obtaining 
the change of variables with the desired properties, moreover when the datum X* is 
limited in regularity. The price to pay is a technical lemma proven in Appendix A. 

Thus we extend the Lagrangian flow Xs(-,t) associated with the solid as a map- 
ping X(-,t) defined in the fluid part. We denote by Y(-,t) the inverse of X(-,t). Then 
we rewrite system (1.1)— (1.9) in a cylindrical domain. For that, another novelty is 
the use of a well-chosen change of unknowns; We introduce the following unknowns 

u(y, t) = R(t) T u{X{y, t), t), p{y, t) = p{X{y, t),t), 

rather than using the whole Jacobian matrix 

u(y,t)=VY(X(y,t),t)u(X(y,t),t), p(y,t) = p(X(y,t),t), (1.15) 

which is done in [12] for instance, or in several papers which only consider a rigid solid 
(see [17], [18], [7] for instance), or simply suggested in [16]. Let us notice that in our 
case the Jacobian matrix VF(A(-, t), t) actually depends on the space variable, and 
thus the use of this classical change of unknowns (1.15) would lead to unappropriate 
complicated calculations and especially it would require more regularity than we ac- 
tually need for the deformation of the solid. 

The corresponding nonlinear system - satisfied by the new unknowns, written in a 
cylindrical domain - is stated in (3.5)-(3.11). The change of variables we have chosen 
enables us to write this system in the simplest form we have found. In particular, the 
equation of velocities (1.4) on dS(t) becomes (3.8) 

3X* 

u = h'+£jf\ X* + — — on aS(0), 
at 

where the datum X* and its time derivative appear in a simple way. The price to pay 
is that we have to study a system in which the divergence of u is not equal to 0. 
The proof of the existence of local-in-time strong solutions is similar to the one done 
in [17]. For proving that the solution so obtained is actually global in time, we show 
that our framework enables us to apply the techniques developed in [7]; In particular, 
we get regularity on the Eulerian velocity w* associated with the Lagrangian mapping 
X*, and we consider an extension of w* to the fluid part. We also quantify the reg- 
ularity needed on this Eulerian velocity, and we observe that the regularity assumed 
on A* is sufficient. 

The choice of the functional framework for the deformation of the solid. 
The mapping X* is chosen such that its time derivative (representing a velocity of 
deformation) lies in L 2 (0, oo; H m (5(0))) n H^O, oo; H 1 (5(0))), where m > 3 is an in- 
teger. The regularity H 3 in space is considered in order to make X* and its extensions 
of class C . Besides, the way we treat the nonhomogeneous divergence condition (in 
the proof of local strong solutions) requires such a regularity, in space as well as in 
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time. Moreover, the estimates we obtain in the proof of global existence (see Propo- 
sition 6.3) require an Eulerian velocity w* whose the regularity - roughly speaking - 
corresponds with the one chosen for the deformation velocity (see Lemma 6.2). 

Let us quote other works which treat of systems coupling the Navier-Stokes equations 
with some deformablc structure, like the mathematical analysis of the interactions 
between a Navier-Stokes fluid and an elastic or viscoclastic structure: [1] , [8] , [4] , [5] , 
[6]. For the fluid-solid system we consider in the present paper, the case of weak 
solutions (in 3 dimensions) has been recently investigated in [14]. Our approach looks 
like a recent work of [2] is which the authors consider an elastic structure whose the 
regularity of its deformation is limited. The interest of considering deformations of 
the solid restricted in regularity lies especially in the perspective of a work where the 
deformation of the solid would be considered as a control function. 

1.3. Plan. 

In section 2 we make precise the functional framework, for the unknowns written in 
time-depending domains and for the datum X* representing the deformation of the 
solid. In section 3 we extend the flow of the solid to the fluid part; It enables us to set 
a change of unknowns and to write the nonlinear system that has to satisfy the new 
unknowns. The linearized system associated with is studied in section 4. In particular 
Proposition 4.4 is used in the next section 5 in order to define a mapping whose a 
fixed point is a strong solution of the nonlinear system. We then prove that for small 
time this mapping is a contraction in a ball chosen large enough. Section 6 is devoted 
to prove the main result, that is to say that the local strong solution is actually global 
if we assume that the data are small enough. Finally, technical lemmas used before 
are proven in Appendixes A and B. 

2. Definitions and notation. 

We denote by 

T = J"(0) and S = 5(0) 

the domains occupied at time t = by the fluid and the solid respectively. We assume 
that S is regular enough. Note that the boundary of T is equal to dO U OS. We set 
for all t > 

S* (t) = X* (S, t) and T(t) =0\ S*(t). 
Let be T e [0,+co]. We set 

S$=Sx(0,T), Q°=J-x(0,T), 

and 

Q T = (J F(t)x{t}. 

tG(0,T) 

In order to deal with some functional spaces, we use the notation 

L 2 (il) = [L 2 (fi)] 3 , H s (0) = [H s (ft)] 3 , 

for all open subset ft C O. Nevertheless this type of notation will be also used for other 
multidimensional spaces (for tensors) like [L 2 (f2)] 3x3 , [L 2 (f2)] 3x3x3 , [L 2 (fi)] 3x3x3x3 , 
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[H s (fi)] 3x3 , [IT(ft)] 3x3x3 or [H s (fi)] 3x3x3x3 , without ambiguity. Let us now make 
precise the functional spaces that we will consider in order to look for strong solutions 
to Problem (1.1)-(1.9). 



2.1. Functional setting for the unknowns. 

The velocity u will be searched in the following functional space 

U(0, T; F(t)) = L 2 (0, T; H 2 (F(t))) fl H x (0, T; L 2 (F(t))) f] C([0, T];H\F(t))), 
that we endow and define with the norm given by 



u 



'u(0,r;^(i)) _ / ll u (')*)llH a (J 7 (i)) ( i* + 



du , . 



dt+ sup |K-,£)|| H i ( w t)) 
L2(^(t)) te[o,T] 



Analogously we can define spaces of type H Sl (0, T; W 2 and IT 1 (0, T; H S2 (fi(i))) 

for all time-depending domain 0(f), where Si and S2 are non-negative integers. The 
pressure p will be searched in L 2 (0, T; H 1 (J r (i))); At each time i it is determined up 
to a constant that we fix such that Jj?u\ P = 0. Thus in particular from the Poincare- 
Wirtinger inequality the pressures P which lie in J- can be estimated in H 1 (J r ) as 
follows 1 

\\P\\^{r) < C\\VP\\v ( ry 

The same estimate will be considered for other functions which play the role of a 
pressure in J-(0). More standardly, we set 

&'HQt) = l2 (0, T; H 2 (^)) n H^O, T; L 2 (J")), 
and we keep in mind that we have by interpolation 

We also set 

H = {(j) G L 2 (C) | div <j) = in O, D(<f>) = in S, • n = on 90} , 
V = {</> G H^O) | div = in O, D(<f>) = in S, <t> ■ n = on 50} . 

2.2. Functional setting for changes of variables. 

For all integer m > 3, we consider deformations of the solid X* which lie in a space 
that we denote by W m (0,T; l S), and that we define as 

X*:Sx[0,T}^ (J S*(t) x {t}, — G L 2 (0 > T;H m (5))nH 1 (0,T;H 1 («S)) 

te[o,T] 

Note that the mapping Ids lies in this space. We take into account the initial condi- 
tions that we assume on X* by considering the displacements X* — Ids in the space 
W^(0,T;S) that we define as follows 



X* - Ids G W^(0,T;S) 



X* G W m (0,T;S), 

X*(y,0) = y, — (y,0) = Vy G S. 



1 In the following the symbols C and C will denote some generic positive constants independent 
of time and the unknowns. 
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Wc endow it with the norm 



\\X* - Ids||wy(o,T;S) : 



dX* 



dt 



L 2 (0,T;H™(S))nH 1 (0,T;H 1 (S)) 



Notice that we have for T < oo 



W m (0,T;S) H 1 (0,T;H ra (5))nH 2 (0,T;H 1 (5)). 



Thus, for more clarity, we set 



VV ro (S£) 

Hm {Qt) 



H^O, T; H Tn (5)) n H 2 (0, T; U\S)), 
H x (0, T; H m (J-)) n H 2 (0, T; H 1 ^)), 
L 2 (0, T; U m (S)) n H x (0, T; tt\S)), 
L 2 (0, T; H m (J")) n H^O, T; H 1 ^)). 



3. The change of variables and the change of unknowns. 

In order to transform the main system in domains which do not depend on time, 
we first extend to the whole domain O the mappings Xs(-,t) and Ys(-,t), initially 
defined respectively on S and S(t). The respective extensions X(-,t) and Y(-,t) then 
obtained define a change of variables which will be used to set a change of unknowns 
for the main system. The aim is to consider new unknowns (u,p,h',uj) which are 
defined in cylindrical domains. 

3.1. The change of variables. 

Let be T > 0. Let h £ H 2 (0,T ;K 3 ) be a vector and R <E H 2 (0, T ; M 9 ) a rotation 
which provides an angular velocity to whose components can be read on 



Since H 1 (0, To; M 3 ) is an algebra, we have to E H 1 (0, To; M 3 ). For a given mapping 
X* e Wo™(0, oo;<S) which satisfies the constraint 



for some T > 0, and which is such that for all t £ [0, T) the function X(-,t) maps 
T onto JF(t), dS onto dS(t), and leaves invariant the boundary dO. For that, let us 
first construct an intermediate mapping. 

Lemma 3.1. Let m > 3 be an integer. Let X* be a mapping such that X* — Ids G 
Wo n (0, oo; S) and which satisfies for all t > the equality 





the purpose of this subsection is to construct a mapping X which satisfies 




in J 7 x (0,T), 
on dS x (0,T), 
on dO x (0,T), 




(3.1) 
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Then for T > small enough, there exists a mapping X G W m (<55 n ) satisfying 
detVX = l mJ"x(0,T), 

X = X* ondSx(0,T), (3.2) 

X = R T (Id-h) ondOx(0,T), 

and the estimate 

\\X - Idjr||w m (Q°) < C (\\X* - Id S || Wm(s o ) + II^-'Uh^O.To* 3 ) + IMIh^O.To;!: 3 )) - 

for some independent positive constant C - which in particular does not depend on T. 
Besides, if X\ and X^ are the solutions of problem (3.2) corresponding to the data 
{X* ,h\,R\) and (X* ,/i2,R.2) respectively, with 

hi(0) = ft 2 (0) = 0, Ri(0) = R 2 (0)=I ffi 3, frl(0) = /4 wi(0)=wa(0), 

then the difference X2 — X\ satisfies 

11^2 - ^l|lwm(0§.) - C (\\h' 2 - /»i||Hi(0,To;R 3 ) + 11^2 - ^1 !|hi(0,T ;M 3 )) j 

where the constant C does not depend on T . 

The proof of this lemma is given in Appendix A. For all t <G [0, T) we denote by 
Y(-,t) the inverse of X(-,t). We now directly set 

X(y, t) = h(t) + R(t)X(y, t), (y, t) e^x (0, T), 
Y(x, t) = Y(R(t) T (x - h(t)), t), x€ T(t), t e (0, T). 

3.2. Rewriting the main system in cylindrical domains. 

Let us transform system (1.1)— (1.9) into a system which deals with non-depending 
time domains. For that we make the change of unknowns 

u(y, t) = R(t) T u(X(y, t), t), u{x, t) = R(t)u(Y(x, t), t), , - 

p(y, t) = p(X(y, t), t), p(x, t) = p(Y(x, t),t), ^ 



for x e J-{t) and y G T . The change of variables X used to define this change of 
unknowns has been constructed in the previous subsection. We also set 

h'(t) = R(t) T h'(t), Q(t) = R(t) T Lo(t). (3.4) 

Remark 1. Let us notice that if h! and u> are given, then by using the second 
equality of (3.4) we see that R satisfies the Cauchy problem 

A (R ) = S (RD)R = RS(a) 







-U) 3 


with §((D) = 


■ 







\—u>2 





R(t = 0) — I R 3, \—u}2 

So R is determined in a unique way. Thus it is obvious to see that in (3.4) the vectors 
h' and lo are also determined in a unique way. Moreover, since we have 

u(x, t) = R(t)u(Y(x, t),t), p(x, t) = p(Y(x, t),t), 
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and since the mapping Y depends only on h, u> and the datum X* , we finally see that 
if (u,p,h',Q) is given, then [u,p,h' ,uS) is determined in a unique way. 

Using the change of unknowns given above by (3.3) and (3.4), system (1.1)— (1.9) 
is rewritten in the cylindrical domain T x (0, T) as follows 

Qu ~ 

— - vLu + M(u, h',u>) + Nu + u)(t) A u + Gp = 0, y G T, t e (0, T), (3.5) 
at 

diY u = g, yeJ, te(0,T), (3.6) 

u = 0, ye dO, t€(0,T), (3.7) 
dX * 

u = h'(t)+u>(t)AX*(y,t) + —-(y,t), ye OS, t e (0,T), (3.8) 

at 



Mh"(t) = - a(u,p)VY(X) T ndT — AIuj(t) A h'(t), t G (0,T), (3.9) 

JdS 



I*(t)Q'(t) = - / X*(y,t) A (a{u,p)VY{X) T n) AT 
JdS v ' 



-r'(t)u)(t)+r(t)Q(t)AQ(t), tG(0,T), (3.10) 

u(y,0) = u (y), y G T , h'(0) = hi G M 3 , w(0) = cj G M 3 , (3.11) 
where [•]$ specifies the i-th component of a vector 

[Lu]i(y, t) = [Vu(y, t)AY(X(y, t), *)]« + V 2 ^(y, t) : (vf VY T ) (X(y, f), t), (3.12) 

Mfafc', *)(!/,*) = -Vu(y,t)VY(X(y,t),t) (h'(t) + w A X(y, t) + ^~{v,t)\ , 

(3.13) 

Nu(y, t) = Vu(y, t)VY{X{y, t),t)u{y, t), (3.14) 
Gp(j/, t) = VY(X(y, t),t) T Vp(y, t), (3.15) 

a(u,p)(y, t)=v (Vu(y, t)VY(X{y, t), t) + VY(X(y, t), t) T Vu(y, if) - p(y, t)l v 
and 

g(y, t) = trace (Vu{y, t) - VY (x(y, t), i))) 

- V«(y,t) : H H 3 - VF (x(y^),t) T ^ . (3.16) 

Notice that we have 

cof(VX) = det(VX)VY(X) T = VY{X) T , 

and then from the Piola identity we can actually express this nonhomogeneous diver- 
gence term as g = div G, where 



G(y, t) = (l R3 - VY(X(y, t), t)) u(y, t). 
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We can verify that G G H 2,1 (Qy), and under the condition (3.7) we have G = on dO. 

Searching for solutions (u,p,h' ,u) to system (1.1)— (1.9) in the space 

W(0,T; F(t)) x L 2 (0,T;H x (J-(i))) x 1^(0, T;R 3 ) x H^O.T;! 3 ). 

is equivalent to searching for solutions (u,p,h',ui) to system (3.5)— (3.11) in 

H 2A (Q° T ) x L 2 (0,T;U\T)) x 1^(0, T;R 3 ) x 1^(0, T;R 3 ). 

The care of verifying the calculations is left to the reader; For that, the main tools 
are the chain rule and change of variables formulas in integrals (see [11] for instance). 

In order to consider a linearized system, we rewrite the nonlinear system (3.5)- 
(3.11) as follows 

-^-uAu + Vp = F(u,p,h',uj), inJ"x(0,T), (3.17) 
at 

div u = div G(u,h',ui), in T x (0,T), (3.18) 

u = 0, iadOx(0,T), (3.19) 
u = h'(t)+Q(t)Ay + W(w), (y, t) G dS x (0, T), (3.20) 

Mh" = -[ a(u,p)ndT + F M (u,p,h',Q), in (0,T), (3.21) 

JdS 

I Q'(t) = - yAa(u,p)ndT + Fi(u,p,h',Q), in (0,T), (3.22) 

JdS 

u(y,0) = u (y), y G T, h'(0) = hi G M 3 , w(0) = w g R 3 , (3.23) 

with 

F(u,p, h',£j) = v(L - A)u. - M(u,h',u)) — Nu — (G - V)p - <D A u, 
G(u, h', u) = (Irs - VY(X(y, t), *)) u, 

W(w) = u A (X* - Id) + 



9i 

F M (u,p, h!,Q) = -MwAh'{t) 



-v J (vu (VY(X) -I R3 ) + (VY(X) -I K 3^) T V5 r j \7Y(X) T ndT 



a(u,p) (vf(X)-I R3 ) T ndr, 
/as v ' 

Fi(u,p, h',u>) = - (I* -I )Q' -!*'& + I*u> A Q 



v J y A (Vu (VY(X) - Irs) + (VY(X) - I R 3) T Vu T ) VY(X) T ndT 
V A (<r(ii ) p)(VF(X) - I R 3) T n"l dr 



+ / (X* - Id) A (a(u,p)S7Y(X) T n) dr. 
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Remark 2. An important remark is the following: Since systems (1.1) —(1.9) 
and (3.17)-(3.23) are equivalent, and since under Hypothesis H2 the compatibility 
condition is satisfied for system (1.1) —(1.9), in system (3.17)-(3.23) the underlying 
compatibility condition enables us to have automatically the following equality 

[ G{u, h', u>) ■ ndT = [ W{oj) ■ ndT 

JdS JdS 

as soon as u = on dO . 

4. The nonhomogeneous linear system. 

Let F, G, W, Fm and Fj be some data. We assume that G satisfies the homogeneous 
condition 

G = on dO 
and also the compatibility condition 

/ G-ndr= / W-ndr. 

JdS JdS 

We assume that 

FeL 2 (0,r;L 2 (J)), GeH 2 - 1 ^), 

W e L 2 (0, T; U 3 / 2 (dS)) n H x (0, T; H 1 / 2 ^)), 

F M e L 2 (0,T;M 3 ),F 7 g L 2 (0,T;]R 3 ). 

In this section we consider for < dist(<S, dO) the following linear system 



— -isAU + VP = F 1 inJ"x(0,T), (4.1) 
at 

div U = div G, in J" x (0, T), (4.2) 

U = 0, ondOx(0,T), (4.3) 
U = h'(t) +w(t) Ay + W, ye OS, t€(0,T), (4.4) 

Mh"(t) = - [ o-(U,P)ndT + W M , t€{0,T), (4.5) 

JdS 

I u'(t) = - I yAo-(U,P)ndr + Fi, t € (0,T), (4.6) 

JdS 

U(y, 0) = u (y), y g JT, h'(0) = hi € R 3 , w(0) = uj g K 3 . (4.7) 
4.1. A lifting method. 



Let us first eliminate the nonhomogeneous divergence condition: By setting 

U = U — G, P = P, H' = h', O = Cj 
we rewrite system (4.1)-(4.7) as 
dU 

— -uAU + VP = F, inJ"x(0,T), 
ot 

div[/ = 0, inJ"x(0,T), 
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17 = 0, ond0x(O,T), 
U = H'(t) + fi(f) A y + W, y € dS, t g (0, T), 

Afff"(t) = - / P)ndr + F M , t g (0, T), 

I Q'(t) = - [ 2/A<r(*7,P)ndr + Fj, ig(0,T), 

U(y, 0) = uo(y), y e T, JJ'(O) = h x g M 3 , O(0) = w G R 3 , 



with 



<9G 

F = F — — + 2vD{G), VF = W-G, 

F M = F M - 2u I D(G)ndT, F I = ¥ I -2v [ y A F>(G)ndr. 

dS JdS 

We now use a lifting method in order to tackle the non-homogeneous Dirichlct condi- 
tion W on dS and then establish an existence result for the linear system (4.1)-(4.7). 
We split this problem into two more simple problems, by setting 

U = V + W, P = Q + 1T, 

where, for all t € (0,T), the couple (w, tt) satisfies 

-vAw(t) + V7r(t) = 0, in F, (4.8) 

div w(t) = 0, in J 7 , (4.9) 

r(t)=W(-,t), ondS, (4.10) 

= 0, on dO, (4.11) 



wi / 

wl 



and where the couple (V,Q) satisfies 

dV 

— - vlSV + VQ = F, inFx(0,T), (4.12) 

div^ = 0, inFx(0,T), (4.13) 

V = 0, on <90 x (0,T), (4.14) 

V = H'(t) + £l(t) A y, ye dS, te(0,T), (4.15) 

MH"{t) = - I a(V,Q)ndr + F M , t € (0,T), (4.16) 

J O'(t) = - / j/Aa(F,Q)ndr + Fr, te(0,T), (4.17) 

7 05 

%0)=« fl (»)-w(»,0),i(€7 1 ff'(0) = hi g R 3 , O(0) =^ elR 3 , (4.18) 

with 

9w 

F = F — — , W = W, 

F M =Fm+ / cr(w, ?r)ndr, Fj = Fj + / y A ct(w, 7r)ndr. 
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To sum up, we have as right-hand-sides: 




(4 
(4 



19) 
20) 



W = W - G, 





<r(w, 7r)ndr, 



21) 



Fj = ¥j - 2u I y A D(G)ndT + 



y A <r(w, 7r)ndr. 



22) 



Jos 



4.1.1. Stokes problem. 

We now look at the problem (4.8)-(4.11). Let us keep in mind that we have the 
compatibility condition 



Let us set a result of existence and uniqueness in H 2 ' 1 ^!^) x L 2 (0, T; H 1 (J 7 )) for this 
nonhomogcncous boundary problem, which is a consequence of a result stated in [9] , 
Theorem 6.1, Chapter IV. 

Proposition 4.1. 

There exists a unique couple (w,7r) G H 2,1 (Qy) x L 2 (0, T; H^J 7 )) solution of system 
(4.8) -(4.11) for almost all t € (0,T). Moreover, there exists a positive constant C 
such that 

II w IIh 2 .i(Q° ) + l|V7r|| L 2( 0j X;L a (^)) < C (|I^IIh 2 .1(Qt) + H^ r llL2(0,T;H 3 / 2 (9S))nH 1 (0,T;H 1 /2(as)) 

The estimate we give in this proposition is not sharp, but it is sufficient for what will 
follow. 

4.1.2. Semigroup approach. 

We solve (4.12)-(4.18) by proceeding in the same way as it is done in [18]. We project 
the unknown V on the space 

H = {4> e L 2 (C) | div (j) = in O, D{4>) = in S, 4> ■ n = on dO} , 

and we consider 

V = {(j) G H x (0) | div 4> = in O, D{<f>) = in S, cf>- n = on dO} . 
Let us recall a lemma stated in [19, page 18]. 
Lemma 4.2. 

For all 4> G % , t/iere exists 1$ e M 3 and cj^ e M sttc/i £/ia£ 



This result allows us to extend V in 5 and then consider the system in the whole 
domain O. Indeed, for V G %, this lemma gives us two vectors H 1 and Vt such that 




4>{y) = l<f> + w </> A y for all y e S. 



V = H' v (t) + Q v (t) Ay — H'(t) + n(t) A y. 
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Let us now define a new inner product on L 2 (0) by setting 

W^k^o) = / ("0 • 0)dy + as / i>(y) ■ 0(y)dy. 



We recall that p$ > is the constant density of the rigid body S. The corresponding 
Euclidean norm is equivalent to the usual one in L 2 (0). If two functions ijj et <j> lie 
in H, then a simple calculation leads us to 



In order to solve (4.12)-(4.18) we use a semigroup approach. We define 

D(A) = {<p£ U\0) | 4>\jr G H 2 (J"), div = in 0, D(</>) = in 5, ■ n = on dO) 

For all F G we set 

and 




^ / L>(y)ndr + \ 2uIq- 1 I y A L>(V)ndr J A y in 5, 
M J as V ./as / 



AV = PAV, 



where P is the orthogonal projection from L 2 (C) onto %. Then we get a unique 
solution (V,Q,H',fl) in H 2,1 (Q^) x L 2 (0, T; H 1 ^)) x H^O.TjR 3 ) x H x (0,T;K 3 ) by 
following the steps of [18]. 

4.2. The main result for the linearized system. 



Definition 4.3. 

A quadruplet (U,P,h',ui) is a solution of the linear problem (4.1)-(4.7) if there exists 
two couples (V,Q) et (w, 7r) such that (w(i), ir(t)) is the solution of system (4.8)-(4.11) 
(given by Proposition 4-1) for all t £ (0,T), such that (V,Q,h',Q) is the solution of 
the problem (4.12)-(4.18) (given by the semi-group approach above), and such that 

(U,P) = (V,Q) + (w,ir) + (G,0). 



Proposition 4.4. 

LetF G L 2 (0,T;L 2 (J-)), F M G L 2 (0,T;R 3 ) ; Fj G L a (0,T;R 3 ), G G H 2 < X (Q°,) and 
W G L 2 (0,T;H 3 / 2 (&S)) n H x (0, T; H x / 2 (dS)) be given. Let G satisfy G = on dO 
and the compatibility condition 

[ G-ndr= / W-ndr. 

J dS JdS 

Assume that < dist(<S, dO) and that uq G H 1 (J r ) with 



div uq = in J-, uq — on dO, Uo(y) — hi + u>q A y on dS. 



16 



S. COURT 



Then system (4.1)-(4.7) admits a unique solution (U,P,h',ui) (in the sense of Defi- 
nition 4.3) in H^^QO) x L^TjH 1 ^)) x H x (0,T;R 3 ) x H x (0,T;R 3 ). 
Moreover there exists a positive constant K such that 

II^IIh^CQO ) + I|V-P||l 2 (0,T;L 2 (JP)) + ||^'||h1(0,T;R3) + |M|hi(0,T;K3) 
< K ^IwoUhMO) + II^IIh^^QO ) + l|Wj| L 2(o,T;H 3 / 2 (a5))nH 1 (O.T;H 1 /2(a l s)) 
+ ||F||l 2 (0,T;L 2 ) + ||1Fm||l 2 (0,T;R3) + ||^/||l 2 (0,T;B 



The constant K depends only on T, Uq, h\, ujq, and is nondecreasing with respect 
to T. 

Proof. Proposition 4.1 provides us a solution (w, n) G H 2a ((5^) x L 2 (0, T; H^J 7 )) 
for the nonhomogeneous Stokes problem (4.8)-(4.11). Let us recall the expressions 
(4.19)— (4.22) of the quantities which appear in some second members of the system 
(4.12)-(4.18): 

„ ^ d& . _ dw 
F = ¥ — — - + i/AG — — , 
at at 

W = W-G, 

F M = F M - 2v / D(G)ndT + / <r(w, 7r)ndr, 

JdS JdS 

F I = ¥ I -2u yAD(G)ndT+ y A c(w, 7r)ndr. 

JdS JdS 

Then the semigroup approach 4.1.2 gives us a unique solution(V, Q, H', fi) for the 
problem (4.12)-(4.18), with 

V G L 2 (0, T; H 2 (J-)) n C* ([0, T]; H 1 ^)) fl H 1 (0, T; L 2 (J")) , 
VQ G L 2 (0, T; L 2 (J")) , ff'eH 1 (0, T; R 3 ) , fie H 1 (0, T; R 3 ) . 

We get then 

(C/,P) = (y,Q) + (w )7 r), 

so by setting (U, P, h', u>) — (U + G, P, H' , f2) we get a solution for the problem 
(4.1)-(4.7) (which is unique). For the wanted estimate, we first write 

II^1Ih 2 >i(Q0,) < ||C/||h 2 -1(QO ) + !|G||h 2 .1(Q5,) 

< II^I!h 2 .i(q?,) + II w IIh 2 -i(q» ) + I|g|| H 2-i(qo ) 

and 

l|VP||L 2 (0,T;L 2 (jr)) < ||VQ|| L 2( 0! T;L 2 (JC-)) + ||V7r|| L 2( 0)T . L 2pr)). 

Then we use the estimate of Proposition 4.1 to get 

II^IIh 2 '!^) + I|vp||l 2 (o,t ; l 2 (^)) < II^IIh 2 ^(q^) + II v< 5IIl 2 (o,t;l 2 (^)) + II^IIhs-mq^) 

+ C (|I^IIh 2 . 1 (Q5,) + l|Wj| L 2( 0jT . H 3/2( aS ) )nH l(0^T;H 1 / 2 (aS)) 

It remains us to use an estimate of the semigroup theory for estimating ||^||h 2 . 1 (q!J,) + 
II VQI|l 2 (o.t ; l 2 (.f)): an d to use again the estimate of Proposition 4.1 to conclude. □ 
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5. Local existence of strong solutions. 
5.1. Statement. 

Theorem 5.1. 

Assume that X* — Ids G W™(0, oo; S) satisfies the hypotheses HI — H4, with m > 3. 
Assume that < dist(<S, <90), and that uq £ H 1 (J r ) satisfies 

div Mo = in T , uo = on dO, uo(y) = hi + luq A y on dS. 

Then there exists Xb > such that problem (1.1)— (1.9) admits a unique strong solution 
{u,p, h, uj) in 

U(0,T o ;F(t)) x L^O^qjH 1 ^))) x H 2 (0,T ;R 3 ) x 1^(0, T ;R 3 ). 

Moreover, if we assume that, for all t £ [0,Tq), dist(iS(i), dO) > r\, then we have the 
alternative 

(a) either Tq = +oo (that is to say the solution is global in time) 

(b) or the function t \— > \\u(t)\\ni(F(t)) * s n °t bounded on [0,Xo). 



Remark 3. The existence of a local strong solution for system (3.5)-(3.11) is 
going to be obtained by a fixed point method for some time Tq small enough. This sys- 
tem is equivalent to system (1.1) -(1.9), up to the change of variables whose existence 
is conditioned by an other time T small enough (see Lemma 3.1 and the change of 
unknowns given in (3.3)J. So, by reducing the existence time T to T, we can get the 
desired local strong solution for system (1.1)— (1.9) . 

5.2. Proof. Let us set 

H T = {(U,P,H',Q.) £ R 2A (Q° T ) x L 2 (0,T;H 1 (J')) x 1^(0, T;]R 3 ) x 1^(0, T;R 3 ) | 
U\ do = in dO x (0,T)} . 

The equivalence of the solutions of systems (1.1)— (1.9) and (3.5)-(3.11) has been 
explained in section 3.2. A solution of system (3.5)-(3.11) is seen as a fixed point of 
the mapping 

(V,Q,K',m) i ^ {U,P,H',Q) 

where (U,P,H' satisfies 
dU 

— - vAU + VP = F(V, Q, K', m), in T x (0, T), 
at 

div U = div G(V, K',zu), in T x (0, T), 

U = 0, mdOx(0,T), 
U = H'{t) + Q{t) A y + W{w), {y, t) £ dS x (0, T), 

MH" = -[ o-(U,P)ndT + F M (V,Q, K' ,w), in (0,T) 

JdS 

I a'(t) = - [ y Aa{U,P)n&T + Fi(V,Q,K\m), in (0,T) 

JdS 
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U(y, 0) = uo(y), ye J 7 , h'(0) = hi G M 3 , w(0) = w e K 3 . 

The expressions of the right-hand-side are given by 

F(V, Q, K', w) = p(L ( k>,*,) - A)V - M (Jf / iW) (V; X', ro) - N ( ^ iCT) y 

-(G (ff , )-V)Q- W AF, (5.1) 

G(V, A", ro ) = fl R3 - Vy(X(y, i), tj) V, (5.2) 



»5 



<9X* 

W(m) = zuA (X* -Id) + — — , (5.3) 

at 

F M (V,Q,K',m) = -Mm A K'(t) 

-v J (WY(X) - Ir:^ + (WY(X) - I M s) T W T J WY(X) T ndT 

-f a(V,Q)hY(X)-I s A T ndr, (5.4) 
Jds v 7 

Ff(v,Q, a',to) = -(/* -7 )Q'-r'o + rr2Afi 

-i^y y A (w (VY(1) - Irs) + (VY(X) - Irs) t W T ) VY(X) T ndT 
y A a(V, Q) (vY(X) -l R3 ) ndT 

i 

+ / (X*-Id)A (t(V,Q)VY(X) T n)dT. (5.5) 

The mapping X is given by Lemma 3.1, with (K', zu, X*) as data. For the expression 
of F(V, Q, K' , w), let us recall that 

[L(jf, ro )(F)]i(y, t) = \VV(y, t)AY(X(y, t),t)} t + V 2 Vi(y, t) : (vfvf T ) (X(y,t),*), 
M (jr / )07) (V; A', a7)(y, t) = -W(y, t)VY{X{y, t),t) (^K'(t) + w A t) + ^(y, t)^ , 

N (K , iVj) V(y,t) = VV(y,t)WY(X{y,t),t)V(y,t), 
G {K , M Q{y,t) = VY(X(y,t),t) T VQ(y,t). 

5.2.1. Preliminary estimates. 

The estimates given in the lemmas below are not necessarily sharp, but they are 
sufficient to prove the desired result. 

Lemma 5.2. 

There exists a positive constant C such that for all (V,Q, K' , w) in Mt we have 



(A L) y|| L 2( 0iT;L 2pr)) < C ||^|| L 2( T .H2(^)) 



X 



vf(x)vy(x) T -i R3 || L oc (0iT: H m - 1( ^) ) + ||Ay(x(.,t),oilL~ ( o,T;H--(^) ) ), 

AY(X(-,t),t)\\ L oa^ 0T . Hm -2^ < C||VF(X) - I R 3||l~(0,T:H™-i(^'))II V ^(^)IIl° c '(0,T;H™-1(J : -)), 

(V - G)Q|| L 2 (0iT . L 2 ( _ F)) < C\\VY(X) - I R 3|| L oo( ^ T;H ™-i(^)) ||V<5|| L 2 (0 T;L 2 ( _ F)) ■ 
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Proof. Given the regularities stated in Lemma 7.4 and the continuous embedding 
H™" 1 ^) ^ L 00 (J r ), the only delicate point that has to be verified is AY(X) e 
L oc (0,T;H m - 2 (.F)). For that, let us consider the i-th component of AY(X); We 
write 



AYi(X(;t),t) = trace (V 2 y(X(-, t), t)) 



with 

V 2 K: 



= (v (vy(x (■,*),*))) vf(A > (-,t),t) 
= (v (vy J (l > (.,t),o - i R3 )) VY(x(-,t),t), 

and we apply Lemma 7.1 with s = m — 2, /Lt = and K = 1 to obtain 

||Ay i (x(. J t),t)|| H m- aw <c||vy(x)-i R 3||H—i ( ^||vy(x)|| H m-i(^. 

□ 

Corollary 5.3. 

There exists a positive constant C such that for all (V, Q, if', ro) in Mt we have 

||(A - L) ^ / || L 2( T;L 2(_ 7 r )) < CVT \\V\\ L 2^ T . n2 ^ X 

(||Vy(l) - ^Hhmo^h—^)) (l + ||Vy(X)|| L oo ( o, T; H™- 1( ^)))) , (5.6) 

j|(V - G)Qj| L 2 (0 T;L 2(.7r)) < CVr\\WY(X) - Ir3||h1(0,T;H™-1(^)) II V QIIl 2 (0,T;L 2 (.F)) ■ 

(5.7) 

Proof. Since VY(X(;0), 0) - I R3 = 0, we have 

||Vy(X) - Ik-'IIl^co.TiH™- 1 ^)) - VT\\VY(X) - I R 3 ||H 1 (o,T ; H m - 1 (j r ))- 
The following quadratic term is treated as follows 

VY(X)VY(X) T - I R 3 = (VY(X) - I R3 ) VY(X) T + (vY(X) T - I R .) . 

□ 

Lemma 5.4. 

There exists a positive constant C such that for all (V, K' , w) in H 2 ' 1 (<5y)xH 1 (0, T; K 3 ) x 
H^TjR 3 ) we have 



\M(V,K',m)\\ L2{0T .^ y) < CT 1 /iO||vy(X)|| L c. (0jT . ir n-i ( ^ ) 

/ llL°°(0,T;H 1 (J r )) II " llLa(0,TjH3(^)) ' 



Ft Y 

K' + wAX + ^- 
ot 



L^O.T.-H 1 ^)) 



^llL=°fO.T:H i m , l H^llL2(-0.T:H 2 m^ ' ( 5 '^) 



II N11/ IIl 2 (0,T;L 2 (^)) ^ C Tl/10 ||Vy(X)|| L oo (0!T . Hm -i ( ^ )) HVllLifo.TjH 1 ^)) H^Hl2(0,T;H 2 (^)) ' 

(5.9) 

A ^I|l 2 (UT;L 2 (.F)) < CVr|jro|| L oo( ,T;B 3 )||^||l~(0 ,T;L 2 (.F))- (5.10) 
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Proof. Let us recall an estimate proved in [18] (Lemma 5.2) which is still true in 
dimension 3; There exists a positive constant C such that for all v, w in H 2,1 (<5^) we 
have 



\(w ■ V)w|| L 5/2 (0iTjL 2 ( ^ )) < C||u;||loo(o,t,hm^))II w 



1/5 



1 4/5 



L~(0,T,H 1 (J !r )) 11 ^lL 2 (0,T,H 2 (.F))' 

(5.11) 

By wlying U- e Sto a te (5.11) with . - Vrt . - - vnl, (it + n A X + f ) , 
combined to the Holder inequality which gives 

\\(w ■ V)u|| L 2(0,T,L=(^)) < T 1/10 ||(u; • V>|| L 5/2 (0 ,T,L2(^)) ! 

we get 



I^^IIl^O.TiL 2 ^)) - 



^fjil/10 



VF (X) ( A" + tu A 1 + ^ 



Il°°(o,T;H 1 (.f)) II^IIl 2 (o,t ; h 2 (.f)) 



L°°(0,T;Hi(J r )) 



We apply Lemma 7.1 on w with s = 1, [i = rn — 2 and k = 0, and then we obtain 
(5.8). For the estimate (5.9), we proceed similarly; We use the inequality (5.11) with 
v = V and w = \7Y(X)V, and we apply Lemma 7.1 on w with s = 1, fi = m — 2 and 
K = 0. For the estimate (5.10), we simply write 

\\zu A V||l2(o,t ; l 2 (:f)) < C|MIl 2 (o,t ; r 3 ) A V^||l«=(o,t ; l 2 (.f)) 

< Cv / r||n7|| L oo(0,T;R3)||^ / ||L~(0,T;L 2 (J r ))- 



Lemma 5.5. 

There exists a positive constant C such that for all (V, K' , zu) in H 2,1 (Q!ji)xH 1 (0, T\ K 3 ) x 
R 1 (0,T;M 3 ) we have 

\\G(V, K',m)\\ h a ( Q tT . l p^ < CVr||^||L 2 (0,T;H 2 (Jf))||VF(X) - I R 3!|h1(0,T;H 2 (^)) 5 

||G(y, A",tz7)|| H i (0T . L 2(jr)) < CVT ^|y|| H i(o,T;L 2 (^))l|VF(X) - I R 3||h 1 (o,t ; h 2 (.f)) 

+ l|y||L 2 (0,T;H 2 (J r ))||VF(X) — lR3||H 2 (0,T;L 2 (J r ))nH 1 (0,T;H 2 (J r )) 



Proof For m > 3, VY(X) lies in H 1 (0, T; H 2 (J")). We apply Lemma 7.1 with 
s = 2, /i = and « = 0, and we get 

\\G(y,K',w)(;t)\\^) <c!|y!| H2( ^||vy(x(-,t),t)-/||H 2 ^), 

||G(V, A', ra7)|| L 2( 0iT .H2(jr)) < C ||y|| L 2 (0,T;H 2 (^)) II VF(X) - I R 3 || L ~(0,T;H 2 (^)) 

< C\/T||V|| L 2 (0 ^ t .h 2 (^)) l|VF(X) - Ir 3 ||h 1 (0,T;H 2 (.F))- 

For proving the regularity H 1 (0, T; L 2 (J 7 )), we first write 
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Notice that we have the embedding 

L 2 (0, T; H" 1 - 1 ^)) n H^O, T; L 2 (J")) ^ L°°(0, T; H™/ 2 " 1 / 2 ^)), 
and thus - by applying Lemma 7.1 with s = 0, fj, = m/2 — 1/2 and /c = 1 - the estimate 



at 



< C \\VY(X) - I R 3 ||l°°(o,T;L°°(.f)) 

L 2 (0,T;L 2 (.F)) \ 



ay 



at 



It' 



L~(0,T;H™/ 2 - 1 /2(jr)) 



L 2 (0,T;L 2 (JF)) 

IjV'llLoofO.^HVJ 7 )) 



□ 

Lemma 5.6. 

There exists a positive constant C such that for all {V, Q, K' , w) in Ht we have 



H W/ ( tn )llL 2 (0,T;H 3 / 2 (aS))nH 1 (0,T;H 1 / 2 (S5)) - C (VT|M|h1(0,T;B! 3 ) + lj \\X* - ld S || w o m (0,oo;5) , 

\\F M {V,Q,K',zu)\\ L2(0T . M3} < 

C ^v / Tj|A''j| L o ( 0jT;R 3)||n7|| L oo(o !T;K 3) + (||y||L 2 (o,T ; H 2 (^)) + ||Q||l 2 (o,t ; h 1 (^))) X 
T\\WY(X) — lR3||Hi(o,Tji, ao (as)) ( ||Vy"(X)|| L ~((o,T)x9S) + 1 



lL 2 (0,T;R 3 ) 



< 



C (7 1 ||^*'||l°°(o,T;K 9 ) IMIh^o^r 3 ) 

+VTj|/*'|| L oo( 0iT . R 9)||n7|| L oo( 0iT . R 3) + v / T||/*||l=°(o,T;R 9 )IMI 

+ \/T (||y|| L 2 (0 , T; H 2 (jr)) + ||<5||l 2 (0,T;H1(J c '))) x 
1 + ||Vy(A > )|| L ^ ((0 ,T)x95)) (||VF(1) - lR3||Hi(0,T i L«OS)) + II VX* 



2 

IL°°(0,T;] 



II^*'IIl° c (o,T;R 9 ) < C ; ll-f* - ^o||l~(o,t ; r 9 ) < C? 1 - 
Proof. For the first estimate, we write (for m > 3) 

ll^( ra7 )l!L 2 (0,T;H 3 / 2 (^)) < C||G7|| L 2( 0jT . H 3)||X* -Id|| LO o( 0) T;H3/ 2 (8S)) + C 

ax 



Idas||H l (o,T;L°°(a.s))) , 
(5.12) 



ax* 



at 



< C%/T||zz7|| L 2 ( - 0iT;R 3) 



a/: 



c 



L 2 (0,T;H 3 (5)) 



ax* 



at 



L 2 (0,T;H 3 / 2 (c9.S)) 

5 

L 2 (0,T;H 3 (<S)) 



and 
dW{w) 



dt 



L 2 (0,T;H 1 / 2 (05)) 



< C||ro'|j L 2( T;R 3)j|X* - Id|| L =(0,T;H 1 / 2 (a5)) 



+C||tJ7|| L 2(o i T;R3) 
< CVT||n7'|| L 2( ,T;R 3 ) 

+Cv / r||n7|| L 2 (0:T;R ; 



dX* 



dt 

dX 



L°°(0,T;H 1 /2(05)) 
L 2 (0,T;H 1 (5)) 



d 2 X* 



dt 2 



L 2 (0,T;H 1 /2(a S )) 



dt 

d 2 X* 



dt 2 



L2(0,T;H!(5)) 



d 2 X* 



dt 2 



L 2 (0,T;H!(5)) 
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There is no particular difficulty for proving the other two estimates, if we refer to the 
respective expressions of Fm and Fj given by (5.4) and (5.5). However, let us detail 
the terms due to the inertia matrices. We have 



I*'(t) 
and thus 



Ps 



dX" 



X* I 



dX* 
~dt 



• X* -X* 



dX" 



(y, t)dy, 



\r'(t)\ 



lL°°(0,T;K 9 ) 



< c 

< c 



dX* 



dt 
dX' 



L 2 (S) 



I^"*(-,*)IIl2(5), 



dt 



\x* 



L°°(0,T;L 2 (<S)) 



|L°°(0,T;L 2 (5))i 



^o|| L oo( O T . K 9) < T 



lL°°(0,T;I 



□ 



5.2.2. The mapping Af is well-defined. From Remark 2 and from the esti- 
mates of the previous subsection (see Corollary 5.3 and Lemmas 5.4, 5.5, 5.6), we can 
first claim that the assumptions of Proposition 4.4 are satisfied. Then from Proposi- 
tion 4.4 the mapping Af is well-defined. Moreover we have the following estimate 

\\Af(V, Q, K', w)\\n T < 4 0) (l + \\ F (V, Q, K', w)\\ L=(0,T;L 2 (^)) + l|G(V,i^',n7)|| H 2,i ( Q0 ) 

+ \\F(V, Q, K' , ro)|| L 2 (0 ,T;L2(^)) + ^T\W\\^(o,t-s?) 

+ \\Fm(V, Q, K', ^)||l 2 (o,t ;R 3) + \\Fi{V, Q, K', ^)|!l 2 (o,t ;R 3)) , 

(5.13) 

where the constant is nondecreasing with respect to T, and depends on the data 



\X* 



Ids\\w^(0,T;S)- 



||wo||hi(^), |/ll|jR3, |wo| 

For R > 0, we set the ball 

B R = {(U,P,H',Cl) e H T | 

II^IIh 3 . 1 ^") + 1 1 1 1 L a (O.rjH 1 (J 7 )) + WH'WhUq^t-R 3 ) + II^IIh^O.TiK 3 ) < -R} 

which is clearly a closed subset of Hy. The rest of this section is devoted to proving 
that for R large enough and T small enough the ball Br is stable by Af, and Af is a 
contraction in Br. 

5.2.3. Stability of the set Br by the mapping Af. 

We are in position to claim that for R large enough and T small enough the ball Br 
is stable by Af. 

Lemma 5.7. 

Let us assume that T < 1 and R > 1. There exists a positive constant C , which does 
not depend on T or R, such that for (V, Q, K 1 , w) G Br we have 

\\F{V,P,K',w)\y {0T . W)) < CT l ' w R\ 
\\G(V,K',m)\\ H2 , 1{Q o T) < CVTR 2 , 
||^m(^P,^',^)|| L 2 (0iT . k3) < CVTi? 3 , 
||iM^#>)|| L2 ( ,T; R 3) < CVtR 3 . 
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Proof. These estimates follow from Corollary 5.3 and Lemmas 5.4, 5.5, 5.6 com- 
bined with the estimates (7.6) and (7.12) (given in Appendix A). □ 

By combining Lemma 5.7 and the estimate (5.13), we have for R large enough 
(R > C^ ] ) and T small enough 

N{B R )dB R . 

5.2.4. Lipschitz stability for the mapping Af. 

Let {Vi.PiJ^wi) and (V 2 , P 2 , K' 2 , w 2 ) be in B R . We set 

(U 1 ,P 1 ,H[,fl 1 ) =Af(V u Q u K' 1 ,m 1 ), {U 2 ,P 2 ,H' 2 ,Q 2 ) = Af(V 2 ,Q 2 ,K' 2 ,zu 2 ), 

and 

U = U 2 -U U P = P 2 -P U H' = H 2 -H{, fi = f2 2 -fii, 

V = V 2 -V U Q = Q 2 -Qi, K' = K 2 -K[, w = w 2 -w x . 

We also denote by X\, VYi(Xi) the mappings provided by Lemma 3.1 with (K[, X*) 
as data, and similarly X 2 , S7Y 2 (X 2 ) the mappings provided by {K' 2 ,w 2 ,X*). 
The quadruplet (U,P,H',&1) satisfies the system 

dU — 

— — uAU + VP = F, inJ"x(0,T), 

div£/ = divG, inJ"x(0,T), 

U = 0, in <90 x (0,T), 
U = H'{t) + fi(t) Ay + W, (y, t) e dS x (0, T), 

MH" = - I <t(U, P)ndT + F M , in (0, T) 

JdS 

I n'(t) = - I y A a{U, P)ndT + Fj, in (0, T) 

JdS 

U(y,0) = 0, in F, H'(0) = e M 3 , 17(0) = 0e R 3 , 

with 

F_ = F(V 2l Q 2l K' 2 , 2 ) - F(Vi, Qi, Jf(, n x ), 
G = G(V 2 , n7 2 ) - G(Vi, X 1; tui), 
W = W(*E7 2 ) - W(roi) = in A (X* - Id s ), 
FM = i i M(F 2 ,Q 2 ,^,w 2 )- J F M (Fi,Qi,if(,w 1 ) ) 
Fj = Fj(V 2) Q 2) ^,ro 2 )-F J (yi,Q 1 ,i<:(,wi). 

In particular, Proposition 4.4 provides for this nonhomogencous linear system the 
following estimate 

||([/,P,#',ft)|| HT <Gf (llFjl L 2 (0,T;R3) + ||G|| H 2,i(Q0) + ||W|| H i( 0)T .H3/3(as)) 

+ 11-^/ IIl 2 (0,T;R 3 ) + II^/||l 2 (0,T;K 3 )) • (5-14) 
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Notice that the right-hand-sides F, G, F m and Fj can be written as polynomial 
differential forms, multiplicative of one of the quantities 

V, Q, K' , w, {Xz-Xi), (v? 2 {X 2 )- Vfi(Xi)). 

For instance, the nonhomogeneous divergence condition G can be written as 

G = (VY 2 (X 2 ) - Vf^XO) V 2 + (Vf^lO - I R3 )V. 

We have in particular 

x 2 (;0)-x 1 (;0) = o, vy 2 (x 2 (-,o) ) o)-vyi(x 1 (- ) o),o) = o. 

The mapping VY" 2 (X 2 ) — VYi(Xi) satisfies the estimate (7.13) stated in Lemma 7.4, 
which is useful in order to make TV a contraction. More specifically, the estimates 
(7.7) and (7.13) are rewritten as 

11^2 - ^lllH^H^nH^H 1 ) < C (||i4r'||Hi(o,To;K3) + IMIh^O.To.-R 3 )) , 
\\VY(X) 2 - VF^X^IlHifH^-^nH^L 2 ) < C (II^''I|h1(0,To;K 3 ) + |MIh1(UT ;R )) • 

Then we state the following result, which can be proven with the same techniques 
that have been used for obtaining Lemma 5.7. 

Lemma 5.8. 

For R large enough and T small enough, there exists a positive constant C - which 
does not depend on T or R - such that 

IFIIlwW)) ^ CT 1/W R 2 \\ (V, Q,K' ,w)\\u T , 

II^IIh^QO) - (|I^IIh2.1(Q0) + \\K'\\-ti^(Q.T;K 3 ) + IMIh^O.TjR 3 )) . 

\\ W \\ni(0,T;V/i(dS)) ~ C^\\ m \\HH0,T-,S. 3 ), 

IF^L(o,T;R 3 ) < CVTR 2 \\(V,Q,K',m)\\ MT , 

,w)|| Ht . 

With regards to the estimate (5.14), we deduce from this lemma that for T small 
enough the mapping N is a contraction in Br. Then the first part of Theorem 5.1 is 
proven. The proof of the alternative announced in a second time is classical; See [7] 
(section 3.3) for instance. 

6. Global existence of strong solutions. 
6.1. Statement. 

Theorem 6.1. 

Assume that the hypotheses in Theorem 5.1 hold true. Assume that ||uo||H 1 (j r ); I^-iIr 3 
and \bj$ | R 3 are small enough, and that the displacement X* — Ids * s small enough in 
W™(0, oo; S), for m > 3. Then we are in the case of the assertion (a) in Theorem 
5.1, that is to say that the strong solution of problem (1.1)— (1.9) is global in time. 
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6.2. A preliminary lemma. 



Lemma 6.2. 

Let X* - Id s e W,f (0,T;«S) suc/i that for all t g [0,T) t/ie mapping X*(-,t) is a 
C 1 -diffeomorphism from S onto S*(t). Then the function defined by 



w*(x*,t) 



dX* 
~dt 



(Y*(x*,t),t), x*eS*(t), te[0,T), 



satisfies 

w* g L 2 (0,T;H 3 ( l S*(t)))nH 1 (0,T;H 1 (5*(i)))- 
Moreover, ||w*||L 2 (o,T;H 3 (S*(t)))nH 1 (0,T;H 1 (S*(i))) is an increasing function of 
dX* 

L2(H3(5))nHi(Hi(S)) 

dX 



Ol- 



id tends to when 



dt 



||VF*(X*)|| LO o (H2 ( S )) 5 ||detVX*(-,t)|| Lo o (L oo (s)) , 
goes to 0. 



L 2 (H 3 (5))nHi(Hi(5)) 

The proof of this lemma is given in Appendix B. The aim of this lemma is to 
show that by assuming smallness on \\X* — Ids\\w m (o.oo-,s)i wc impose automati- 
cally smallness on the velocity w* in L 2 (0, T; H 3 (5*(t))) n rl 1 (0, T; H 1 (S* (t))) . Thus 
in the proof of Theorem 6.1 it is sufficient to consider that w* is small enough in 
L 2 (0,T;H 3 (S*(i))) n H x (0, T; H 1 ^*^))) for all T > 0. 

6.3. Sketch of the proof. 

Let us think by absurd. Assume that To < oo. Let us then show that the functions 

t ^ ||u(t)|| H i(^(t)), t^\h'{t)\, t^\u(t)\ 
are bounded in [0,Xb). For that, let us give a first estimate: 

Proposition 6.3. Let (u,p,h,ui) be a strong solution of the system (1.1)— (1.9) 
defined on [0,T ) with T Q > 0. Furthermore assume that there exists r\ > such that 
for all t g [0, T ) 

dist(5(t),aO) > rj. 

Then there exists a positive constant K ( depending on Tq and rj ) such that 

IMk°°(0,To;L 2 (.F(t))) + IMIl^lUTo;!! 1 ^*))) + II h' || L =° (0,T ;K3) + IMIl°°(0,T ;R 3 ) < KC Q , 

C* := exp (i^||w*||L2 (0 ,T ;H3(s.(t)))) x 



(|MIl»(.F) + l /l l| 2 + l^ol 2 + IKIlHi(0,T ;Hi(.S*(t))) (l 



,*l|2 



1/2 



IL2(0,T ;H3(S* (t))) J 



Proof. We need to define an extension to J-(t) of the velocity w(-,t) which is define 
on S(t). For that, let us first define an extension to T(t) of the velocity w*(-,t), defined 
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on S*(t) and whose the regularity has been given in Lemma 6.2; This extension is 
denoted w*(-,t) and is chosen as solution of the following divergence problem 



div w* = 
w*(x*,t) = 
W*(x*,t) = w*(x*,t) 



mm 6 \S*(t), te (0,T ), 

if dist(x*,5*(t)) > ?7 > 0, te(0,T o ), 

if x* G<S*(t), tG (0,T ). 



(6.1) 



A solution of this problem can be obtained by using some results of [9] for instance: 
The nonhomogeneous Dirichlet condition can be lifted (see Theorem 3.4, Chapter II) 
and the resolution made by using Exercise 3.4 and Theorem 3.2 of Chapter III. Then 
this extension w* of the datum w* obeys the following estimates: 



\w*{-,t)\\ H3{j r {t)) < C||w*(-,f)|| H 5/2 (as , (t)) < (5||u;*(-,i)||HS(S'(i)). 



dw* 



dt 



< C 



Hi(^(t)) 



dw* 



dt 



M) 



< C 



dw* 



dt 



(6.2) 
(6.3) 



Hi(S*(t)) 



The constant C does not depend on time, since we have assumed that dist(6>(t), dO) > 
i] > for all t G [0, T ). Then we set as an extension of w in F(t): 

w(x,t) = ~R(t)w* (R{t) T {x - h(t)),t), x G F{t). 

This relation yields the following properties 

div w = 
w = 

w = w 



in F(t), te (0,T ), 
on dO, t G (0,T ), 
on dS(t), t G (0,T ), 



and the following estimates, for some positive constant C independent of time 



MlL 2 (0.T;L 2 (.F(t))) 



W 



lL 2 (0,T;L 2 (.F(t)))' 



l|Vw|| L 2 (^(t)) < C||Vw*|| L2(jE - (t)) , 

(W- V)W|| L 2 (J r (t)) < C||uJ*|| H3( ^ (t)) ||uJ*|| H1(/ . (t)) , 

INIw 1 .-^)) < C W W *\\n3(f(t))' 



dw 



di 



< c 



dw* 



dt 



L 2 (^(t)) 



Ihh^C*)) 



(\h'\ 



(6.4) 
(6.5) 
(6.6) 

(6.7) 



Let us now set v — u — w. The function v satisfies the following system 
Ov cfw 

— + (u-\7)v-vAu + Vp = -(v-V)w-(w-\7)w-—, x G Fit), t G (0,T), 
at at 

(6.8) 

div m = 0, x G J"(t), i G (0, T), (6.9) 



u = 0, xedO, te (0,T), 
u = h'(t) + uj(t) A(x — h(t)), x G dS(t), t G (0, T), 



Mh"(t) = - [ <j{u,p)ndT, t G (0,T), 

J95(t) 



(a - /i(t)) A cr(u, p)ndT, t€(0,T), 



(6.10) 
(6.11) 

(6.12) 
(6.13) 



dS(t) 
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h(0) = h eR 3 , h'(0) = hi G K 3 , w(0) = w g K 3 , 
0) = 110(2;) := uo(^) — w(x, 0), if/. 



(6.14) 
(6.15) 



By taking the inner product of the first equation with v and by integrating on F(t) 
we get after some calculations 
Id 
2ch 



v\hdx + 2vj \D(v)g 3 dx+^-± (\h'(t)\* a ) + ~ (|(^) (t) 

D(«7) : L»(w)dx 

/ ((to ■ V)to) ■ v dx — / — — • v dx. 

JF{t) J Fit) & 



—J w • w — 2f 
2 



.F(t) 

((u ■ V)TtJ) • u dx 



where we use the notation uj = R T w. It follows that there exists C > such that 



2di 

< C 



\ v \k*(F(t)) + 2v \\D(v)\\- L 2 {J r {t)) + 



2 dt 



1 d 

2d7 



y/lu) (t) 





dw 




dt 



L 2 (^(t)) 



+ || (to- V)to|| L 2pr (t)) 



I-CHIIls^*)) 



+ll«llL(^( t) ) (1 + livw||^ (()) ) + | w | 2 (i + |r'||o) ■ 

Using the estimates (6.4)-(6.7) combined to (6.2)-(6.3), we get 



2dt 

< C 



Vlu) (t) 



dw* 



i)t 



,* 11 2 



Hi(5*(t)) 



!H3(5*(t)) 



Klines- 



(S*(t)) 



w *llHi(5*(t)) 



+ (IK'IIL(^)) + Mi? + |/»'||a) (1 + lk*llH3 ( s.(t)) + |r'|| 9 )) 



(6.16) 



Besides, we can extend the velocity field v into <S(£) by setting v(x, t) = h'(t) +w(t) A 
(x — h(t)) for x G <S(i). Thus we have i; G Hq(0) with div v = 0, and the following 
formula 

Vji : V« - 2D(u) : D{v) = -div ((u • V)v - (div i>)v) - (div v) 2 

combined to the Poincare inequality enables us to write 

IMIm(.F(t)) < ||u||hi(0) 

< C||Vu|| L2(0) = 2C||£>(«)|| La(0) - 2C|p( U )|| L2( ^ (i)) . 

Then we can conclude by using inequality (5.12) and the Gronwall's lemma on (6.16). 
□ 

Proposition 6.3 is the analogous adaptation of Proposition 1 and Lemma 4.1 of [16] 
and [7] respectively. The difference with [16] is that in dimension 2 we do not have to 
assume smallness on the data, whereas in our case we need to quantify the regularity 
of the deformation for which we need to assume smallness in dimension 3. In [7], only 
the rigid case is treated (in dimension 2 and 3). From there the rest of the proof is 
quite standard, and we can follow straightforwardly the steps detailed in section 4.3 
of [7] to complete the proof. 
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7. Appendix A: The changes of variables. 

7.1. Preliminary results. Let us recall a result stated in the Appendix B of 
[10] (Proposition B.l), which treats of Sobolev regularities for products of functions, 
and that we state in the particular case of dimension 3 as: 

Lemma 7.1. 

Let s, [i and k in R. If f € H s+A '(J r ) and g € H S+K (,F), then there exists a positive 
constant C such that 

II/S , Hh=(.F) < C\\f\\ ¥L s+^ [:F) \\g\\ Y{s+K( jr ) , 

(i) when s + [i + k > 3/2, 

(ii) with [i > 0, k > 0, 2s + /i + k > 0, 

fmj except that s + fx + k > 3/2 i/ equality holds somewhere in (ii). 



A consequence of this Lemma is the following result. 
Lemma 7.2. 

Let be T > 0. Lei X 6e m Wm(<9!r)j wit/i m > 3. XTien 

cofVX G H 1 (0,T;H m - 1 (J'))nH 2 (0,T;L 2 (J : -)), 
and there exists a positive constant C such that 

||cofVX - iRsllH^H^-^nH^L 2 ) 



< C\\VX - lM 3 \\n 1 (H m - 1 )nlP(V t ) (1 + ||VJf - lM3|| H i(H">- 1 )nH 2 (L 2 

< cVT 



dX 


li + VT 


dX 










dt 


dt 


n m (Q° T )j 




H m (Q° T ) V 







Moreover, if X u X 2 € W m {Q%), then 
||cofVA 2 - cofVA\||Hi(H"«-i)nH 2 (L 2 ) 

< C||VX 2 - VAi|| H i(H™- 1 )nH 2 (L 2 ) ( 1 + IIVXlH + llV^lln^H^-^nH 2 ^ 2 



< cVf 



d{X 2 -X{) 


Ii + Vt 


dX x 


dt 


dt 




n m (Q° T ) \ 





Vt 



H m (Q° T ) 



(7.1) 



(7.2) 



dXi 




dt 


n m (Q° T )j 







(7.3) 



Proof. For proving (7.1), it is sufficient to show that the space H 1 (0, T; H m_1 (J r ))n 
H 2 (0, T; L 2 (J")) is stable by product. For that, let us consider two functions / and g 
which lie in this space. We write 

d(fg) _ df dg 
dt dt 9 1 dt' 
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Applying Lemma 7.1 with s = m — 1 and // = «; = 0, we get 



0f 



< C 



L 2 (0,T;H m - 1 (.F)) 



Of 



at 



IS , IIl~(0,T;H">- 1 (.F)) 



l/ll 



L 2 (0,T;H"'- 1 (^-)) 



L"=(0,T;H~- 1 (J : )) 



and thus /.g e H 1 (0,T;H m - 1 (.F)). For the regularity of fg in H 2 (0, T; L 2 (J r )), we 
consider 

at 2 dt 2 9 J at 2 at at 



with 



|6L 2 (0,T;L-(J)), 



/, g G L°°(0, T; L 00 ^ 7 )), 



||gL°°(0,T;L 2 (.F)), 



because of the embedding H m 1 (J r ) ^ L°°(J r ), so that we get 

d 2 (fg) 



at 2 



G L 2 (0,T;L 2 (.F)) 



and the desired regularity. This shows in particular that H^O, T; H" i - 1 (J r ))nH 2 (0, T; L 2 (J")) 
is an algebra, and we can show the estimate (7.2) by noticing that the cofactor matrix 
is made of quadratic terms (in dimension 3), so that 

||cofVX - I R3 ||H 1 (H m - 1 )nH 2 (L 2 ) 

< C\\VX - lK3|| H i(Hm- 1 )nH 2 (L 2 ) (||VX|| H i( H m-i) nH 2 (L 2 ) + l) i 

< C||VX - lR3||Hi(H m -!)nH 2 (L 2 ) f||VX - lR3||H 1 (H"'- 1 )nH 2 (L 2 ) + ij • 

The argument for proving (7.3) is the same. □ 

7.2. Existence of a change of variables. 

Let be T > T > 0. Let h G H 2 (0, T ; R 3 ) be a vector and R G H 2 (0, T ; R 9 ) a rotation 
which provides the angular velocity uj G H 1 (0,Tq;]R 3 ) given by 



dR 

§ (uj) = — R T , with S(w) 
v ; dt v ' 



-U) 3 uj 2 
u; 3 — (J\ 

—OJ-2 U>1 



We assume that ho = 0, R(0) = Ir3 and we still use the notation 

h'(t) = R{t) T h{t), u{t) = R(t) T uj{t). 
Let us recall and prove Lemma 3.1: 
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Lemma 7.3. 

Let m > 3 be an integer. Let X* be a mapping which lies in Wq™(0, oo; S) and satisfies 
for all t > the equality 



dX* 



(cofVX*)ndr = 0. 



ids 9t 

Then for T > small enough, there exists a mapping X £ W m (Qy) satisfying 



(7.4) 



detVX = 1 

X = X* 

X = R T (Id - h) 



inTx (0,T), 
on dS x (0,T), 
on dO x (0,T), 



(7.5) 



and the estimate 



\X - Id^|| Wm(Q o ) < C [1+ - Ids|| Wm(/ go ) + ||/l'||Hi(0,T ;R3) + ||w|| H i(0,To;R3)J , 

(7.6) 



for some independent positive constant C - which in particular does not depend on 
T. Besides, if X\ and X2 are the solutions - for T small enough - of problem (3.2) 
corresponding to the data (X*,hi,TLi) and (X*, /12, R2) respectively, with 

/n(0) = ha(0) = 0, Ri(0) = R 2 (0) =I R3 , h[(0) = h' 2 (0), wi(0) = wa(0), 

then the difference Xi — X\ satisfies 

11^2 ~ X l\\w m (Q° T ) < <5 (ll^2 ^ ^'iIIh^O.To^ 3 ) + 11^2 - Wi|| H i( 0) T ;R3) J , (7-7) 

where the constant C does not depend on T . 

Remark 4. A mapping X given by the lemma above satisfies X(-,0) = Id jr. 

Proof. Given the initial data X*(-,0) = Ids, K = 0, R(0) = I R3 , h'(0) = hi and 
w(0) = loq, let us consider the system (3.2) derived in time, as 

dVX 



cofVX 



dX dX* 



&t 



= 



at dt 

dX 

— {y,t) = -h'(t)-uj(t)AX(y,t) 
I X(-,0)=Id^. 



in J 7 x (0,T), 
on 55 x (0,T), 
(|/,t)G9Ox(0,T), 



This system can be seen as a modified nonlinear divergence problem, that we state as 

in T x (0,T), 
on 55 x (0,T), 
fe,t)eaOx(0,T), 



div f = /(I) 

ax ox* 



at 
ax 



at 



— (y,t) = -h'(t)-u>(t)AX(y,t) 

{ X(;0)=Idr, 
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with 



f(X) = (l R3 - cofVl) 



<9VA 
dt ' 



If we search solutions to this system which are continuous in space, let us notice (by 
using the Piola identity) that the compatibility condition for this divergence system 
is nothing else than the equality (7.4). 

A solution of this system can be seen as a fixed point of the mapping 



X: W m (Q°) 



W m (Q Q T ) 
X 2 , 



(7.8) 



where X 2 is a solution of the classical divergence problem 



dX 2 _ 8X* 

~df 

dX 2 

dt 



Of 



-h' - wAli 



in J 7 x (0,T), 
on dS x (0,T), 
on dO x (0,T), 



by adding the initial condition X 2 (-,0) = Id jr. A solution of this problem can be 
obtained by using some results of [9] for instance: The nonhomogencous Dirichlct 
condition can be lifted (see Theorem 3.4, Chapter II) and the resolution made by 
using Exercise 3.4 and Theorem 3.2 of Chapter III. Then the solution chosen is the 
one which satisfies the estimates 



dX 2 



dt 



<C[\\f(X 1 )\\v( Hm - H r)) + 



L 2 (H™(^)) 



dX* 



dt 



L 2 (H ™-i/2(as)) 



+ II^'IIl 2 (0,T;R 3 ) + ll^llL 2 (0,T;R 3 )l|A"l|| L oo( Hm -i/2 (0cl ) ) 

dX 



< C II/^OIIl^h-m^)) 



L 2 (H™(S)) 



+ ( II^'I!l 2 (0,T;R 3 ) + ||^||l 2 (0,T;R 3 ) ) 1 + vT 



dXi 
~df 



L 2 (H™(.F)), 



(7.9) 
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and 



d 2 X 2 

dt 2 



<C ||/(Xi)||hi(l»(jo) 



L 2 (Hi(^)) 



ax* 

~dT 



H 1 (H 1 /2(05)) 



+ ||' 1 'I|h 1 (0,T;K 3 ) + IM|H 1 (0,T;R 3 )||^l||L~(H 1 /2(ae>)) + ||w||l°°(0,T;I 

(d X* 
\\f(X 1 )\\ nH v (r)) + 



dXi 




at 






L 2 (Hi/2(SO))/ 



dt 



H 1 (H 1 (S)) 



+ |I^'||h 1 (0,T;R 3 ) + ||w||hi(0,T;I 







dXt 




Vt 




at 



L 2 (H!(^)) 



6X X 




dt 


L~(Hi(^))/j 







(7.10) 



Indeed, let us verify that for X £ W m (Q^), satisfying X(-, 0) = Idjr, we have f(X) € 
L 2 (0,T;H m - 1 (J)) n H'fO.T;^^)). For that, we recall from the previous lemma 
that cofVX G R 1 (0,T;'H m - 1 (J r )) n H 2 (0, T; L 2 (J")), and we first use the result of 
Lemma 7.1 with s = m — 1 and /i = ft = to get 



II/P0IIl 2 (H"'-i(.f)) < C||Ir 3 - cofVX|| L oc( H m-i(^) 



dVX 



< CVT\\l R3 - cofVX|| H i(H—i(^)) 



dt 
dX 



L 2 (Hm -l(^)) 



dt 



L 2 (H™(.F)) 

For the regularity in H 1 (0, T; L 2 (_F)), let us first notice that we have by interpolation 

L 2 (0, T; r-'fj)) n H^O, T; L 2 (J-)) L°°(0, T; H m / 2 -^(J)). 

Then we use Lemma 7.1 with s = and fj, = K = to/2 — 1/2, and the continuous 
embedding H m_1 (J r ) <^-> L°°(J r ) in order to get 

d(f(X)) 



dt 



< C 



L 2 (L 2 ) 

dcofVX 



< CVT 



dt 

L oo( H m/2-V2) 

dcofVX 



<9VX 



9* 



L oc( Hm /2-l/2) 



L 2( H m/2-l/2) 

L«(H«/>-i/i) 



cofVX 



Ir3 — cof VX 





<9 2 vx 




L°°(H™- 1 ) 


<% 2 


L 2 (L 2 )j 







H 1 (H" 





a 2 vx 






at 2 


L 2 (L 2 )j 









Thus, we finally have 

H/(^')llL 2 (H m - 1 )nH 1 (L 2 ) < Cv^HIh 3 - CofVX|j H i(H™-i)nH 2 (L 2 (.F)) 



dX 



dt 



■H m (Q° T ) 

(7.11) 
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The estimates (7.9) and (7.10) combined to (7.11) and (7.2) show that the mapping 
X is well-defined. 

Moreover, the set defined for some R > by 



03fl = {X G W m (Q° T ) | l(-,0) = Idjr, 



9X 



<9t 



<R 



n m (Q° T ) 



is stable by T, for T small enough and R large enough. Notice that 03 # is a closed 

subset of W m (Qy). Let us verify that 1 is a contraction in 03 _r. 

For Xi and X 2 in 03 _r, the difference Z = 1(X 2 ) — 'Z(Xi) satisfies the divergence 

system 



d7 

div — = f(X 2 )-f(X 1 ) 

dZ 

~dt 
dZ 

~dt 



and thus the estimate 
dZ 



in J 7 x (0,T), 
on <9S x (0,T), 

on SO x (0,T), 

< C\\f{X 2 ) - J(^l)|!L2(H'"-i(^))nHi(L 2 (^))- 



L2(H-(^))nHi(HH^)) 



We write 
f(X 2 ) - fix,) 



cofVX 2 - cofVXi 



dVX 2 



Ir3 — cofVXi 



aV^a-Xi) 
St 



By reconsidering the steps of the proofs of the estimate (7.11), and by using (7.3), we 
can verify that for T small enough the mapping % is a contraction in 03 r. Thus T 
admits a unique fixed point in 03 /j. 

For the estimate (3.3), let us just notice that the difference X of two mappings X\ 
and X 2 of 03 - corresponding to the data (X* , hi, Ri) and (X*, h 2 , H 2 ) respectively 
- satisfies the system 

dX 

div — = f(X 2 ) - f(Xi) mTx(0,T), 
BX 

— =0 on<9Sx(0,T), 
dX 

— = -(h 2 - h'i) -{ui 2 -ui 1 )AXi-Q 2 AX on dO x (0, T). 
dt 

Then we proceed as previously, and the end of the proof for the announced estimate 
is left to the reader. □ 



7.3. Lipschitz estimates. 



Lemma 7.4. 

Let be Tq, and T > small enough to define X G W m (Qy) solution of problem (7.5), 
for X* G W m (S^), h G H 2 (0,T ;R 3 ) and R G H 2 (0, T ; M 9 ). Let us denote by Y(-,t) 
the inverse of the mapping X (-,t) - for all t G [0, T). Then we have 

\\VY(X) - lR3|| H l( T . Hm -l(_ 7 r)) nH 2( 0jT . L 2(jr )) 

< C f 1 + ||X* - Id s \\ Wm ( S o + II/i'IIhvo,^^ 3 ) + \\^\\n 1 (a,T X->) ) • (7-12) 
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Let Xi, X 2 € WmiQ^r) ^ e ^ e solutions of problem (7.5) - for T small enough - with 
(X*,ft, 1 ,R 1 ) and (X*,/i 2j R-2) a s data respectively. Then, if we denote byYi(-,t) and 
Y2{-,i) the inverses of X\(-,t) and X2(-,t) respectively, we have 

||VF 2 pT2) — Vyi(Xi)|| H i(o,T;H TO - 1 (^"))nH 2 (0,T;L 2 (^")) 

<C(\\h' 2 — ^illH^O.ToiM 3 ) + 11^2 - ^iIIhHO.Tq^ 9 ) ) • ( 7 -13) 



Proof. By considering the first equality of Problem (7.5), we have the equality 



VY(X(-,t),t) = ^J^;^ = cofVX(-,tf , 
detVA (■, t) 

and so Lemma 7.2 combined to the estimates of Lemma 7.3 can be applied. □ 

8. Appendix B : Proof of Lemma 6.2. 

Let us use a result given in the Appendix of [3] (Lemma A. 4), which treats of regularity 
in Sobolev spaces for composition of functions: There exists a positive constant C such 
that for all t e (0, T) we have 



lk*(-,t)l|H3(S.( t) ) <c 



< c 



dX* 



dt 

8X' 



dt 



(■,*) 



r*(^)iiH3 (s . W ) + i 

hs C s) inf |detVY*(z*,t)| 1/2 



H3(S) 



(\\Y*(;t)fn*(s* { t)) + l) \\dcWX*(-,t)\\^ {sy 



Let us notice that by using the change of variables induced by X*(-,i), we have 



|y*('^)IIl 2 (s*(*)) 



\Y*(x*,t)\l 3 dx* 

S*(t) 

y\ 2 deWX*(y,t)dy, 



l|vr*(-,t)!| L2(5 . (t)) < !|vy*(A*(y,t),i)|| L3(5) ||dctvA*(.,t)||^ 2 (5) . 

The following equality 

V 2 y*(X*(-,£),t) = (V(VY*(X*(-,t),t)))VY*(X*{-,t),t) (8.1) 

yields 

l|v 2 r*(-,i)|| L2(5 . (t)) <c||vy*(x*(- ) i),t)llHi(s)l|vy*(x*(-,i),t)llL-.(s)l|detvx*(-,t)| 

Moreover, by applying Lemma 7.1 with s = 1, /i = and « = 1, the equality (8.1) 
implies 

j|v 2 y*(x*M),t)j| H1 ( 5 ) < q|vy*(x*(.,t),t)||^ cs) . 

The following equality 

v 3 Y*(x*(-,t),t) = (v 2 (vr(r(.,t),t))) (vY*(x*(-,t),t)f 
+ (v (vr(r(.,t), *))) (v 2 y*(A*(-, t), *)) 



1/2 
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combined with the previous estimate enables us to obtain 



l|v 3 y*M)llL 2 ( 5 . W ) <q|vr*(x*(.,i),i)||^ (<s) ||dctvx*(.,t)||^ 2 (5) . 



Finally we get 

l|w*|| L 2(H3(5*(t))) < C 



dX* 



at 



|detVX*(-,t)| 



L=(H3(^)) 



1/2 

L°°(L°°(S)) 
3 



3/2N 



1+ ||detvx*(-,t)|| L ^ (L ^ (5)) ^||vy*(x*)||2i (H2(5)) 



fc=l 



For the regularity of w* in H 1 (0, T; H 1 (5*(t))), we estimate 



(•,*) 



9* 

and we calculate 



<q|detVX*(-,i) ll^s) 



L 2 (5*W) 



9w* 



(9/ 



(**(•,*), t) 



L 2 (S) 



— (X*(y,t),t) = -^j-fo.t) - -^-(| /) i)Vr(r(|/ 1 t) 1 t) — (i/,*), y e 5. 



Thus we have 
9w* 



at 



<C||detVX*(.,t)||^ 2 (Loo(5)) x 



L3(L*(S»(t))) 



a 2 X* 



+ 



avx* 



at 



L 2 (L 2 (5)) 



at 2 

vy*(x*)|| L ^ ( H 2(< s)) 



L=(L=(S)) 



ax* 




at 


L~(H2(S))j 



Finally we control in L 2 (0, T; H x (5*(t))) by writing 

f)X7w* f) 2 X7 X* r)X7 2 X* BY* 

-^-(*V) = ?-^(Y*(x*,t),t)VY*(x*,t) + ^ r (y*( i r* ) t) ) t)Vr( 5t * 1 t)^-(i* ) () 

avx* av\^* 
+ -^-(^*(a:*,*) ) t)-^-(a:*,t), G $*(*), 

-^(X%y,t),t) = ^^(y 7 t),t)\7Y*(X*(y,t),t)-^^(y,t)(\7Y*(X*(y,t) 7 t)r^ r (y,t) 



dVX* 



Of 



a 



+—~r-M ~r (vy*(x*( y ,t),t)) - v(vr*(x*M),t)) , </e5 



at 



and by estimating 



at 



<C||detVX*(.,t)||^ 2 (LOO(s)) x 



L2(L=(S*(*))) 



at 2 



av 2 x* 



at 
avx* 



L 2 (L 2 (S)) 



l|vr*(x*)||^ (H2(s)) 



L 2 (L 2 (S)) 

dX* 



|vr*(x*)|| L o 0(H2 ( 5 ) ) + 



at 



L°°(*P(S)) 



at 



L~(L=(S)) 



|VF*(X*)|| H i(H 2 (s)) 
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